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Bilayer quantum Hall systems, realized either in two separated wells or in the lowest two sub-bands of a wide 
quantum well, provide an experimentally realizable way to tune between competing quantum orders at the same 
filling fraction. Using newly developed density matrix renormalization group techniques combined with exact 
diagonalization, we return to the problem of quantum Hall bilayers at filling v = 1/3 + 1/3. We first consider 
the Coulomb interaction at bilayer separation d, bilayer tunneling energy Asas, and individual layer width w, 
where we find a phase diagram which includes three competing Abelian phases: a bilayer-Laughlin phase (two 
nearly decoupled iz = 1/3 layers); a bilayer-spin singlet phase; and a bilayer-symmetric phase. We also study 
the order of the transitions between these phases. A variety of non-Abelian phases have also been proposed for 
these systems. While absent in the simplest phase diagram, by slightly modifying the interlayer repulsion we 
find a robust non-Abelian phase which we identify as the “interlayer-Pfaffian” phase. In addition to non-Abelian 
statistics similar to the Moore-Read state, it exhibits a novel form of bilayer-spin charge separation. Our results 
suggest that ;z = l/3 + l/3 systems merit further experimental study. 
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I. INTRODUCTION 

The remarkable experimental discovery of quantized resis¬ 
tance of a two-dimensional electron gas (2DEG) in strong 


perpendicular magnetic fields [1] has revealed many topolog¬ 
ically ordered phases that form due to strong Coulomb inter¬ 
actions in a partially filled Landau level [2]. Some exam¬ 
ples include the “odd-denominator” fractional quantum Hall 
(EQH) states that belong to the sequence of Laughlin [3], 
hierarchy [4, 5] and “composite fermion” [6] states. One 
of their prominent features is the presence of quasiparticles 
(“anyons”) that carry fractional charges [3] and obey frac¬ 
tional statistics [7, 8]. More intriguing, “non-Abelian” quasi¬ 
particles have been proposed to occur in several experimen¬ 
tally observed EQH states in the first excited Landau level. 
Most notably, this is the case with an even-denominator fill¬ 
ing factor zz = 5/2 state [9], believed to be described by the 
Moore-Read Pfaffian state [10-12] that contains non-Abelian 
anyons of the Majorana type [13-15]. 

The aforementioned hierarchies of Abelian and non- 
Abelian states are a priori relevant when the EQH system can 
be described as a single partially occupied Landau level, that 
is, the electrons carry no internal degree of freedom. However, 
“multicomponent” EQH states are ubiquitous; most obviously 
electrons carry spin. While the Coulomb energy scales as 
e^/e^B[K] Ri SOi/Sfr], assuming free electron values for 
the mass and g factor in GaAs, the Zeeman splitting is only 
Ez\ii] ~ 0.3i3[T], suggesting that in many circumstances the 
ground state of the system may not be fully spin-polarized. 
Several classes of unpolarized EQH states have been formu¬ 
lated, including the so-called Halperin [mmn) states [16] and 
spin unpolarized composite fermion states [17-20]. In mate¬ 
rials such as AlAs or graphene, ordinary electron spin may 
furthermore combine with valley degrees of freedom, which 
can change the sequence of the observed integer and EQH 
states [21-30]. 

Here we study an important class of multicomponent EQH 
systems where the internal degrees of freedom correspond to 
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a subband or layer index, generally referred to as pseudo-spin. 
For example, if a 2DEG is confined by an infinite square well 
in the perpendicular z-direction, the effective Hilbert space 
may be restricted to several low-lying subbands of the quan¬ 
tum well (QW). In the most common case, the relevant sub¬ 
bands are the lowest symmetric and antisymmetric subbands 
of the infinite square well that play the role of an effective 
SU(2) degree of freedom. Furthermore, it is possible to fab¬ 
ricate samples that consist of two quantum wells separated by 
a thin insulating barrier. We refer to the latter type of device 
as the quantum Hall bilayer (QHB). The interest in bilayers 
and quantum wells comes from their experimental flexibility 
that allows one to tune the parameters in the Hamiltonian to 
a larger degree than it is possible with ordinary spin. For ex¬ 
ample, in a QHB with finite interlayer distance, the effective 
Coulomb interaction is not SU(2) symmetric. Therefore, the 
“intralayer” Coulomb interaction (the potential between elec¬ 
trons in the same layer) is somewhat stronger than the “inter¬ 
layer” Coulomb (i.e., the potential between electrons in oppo¬ 
site layers). The ratio between the two interaction strengths 
is given by the parameter d/is, the physical distance be¬ 
tween layers in units of magnetic length, which in experiment 
can be continuously tuned. The tunneling energy between 
the two layers (in units of the Coulomb interaction energy), 
can also be tuned. The tunability of interactions 
in quantum Hall bilayers and quantum wells can give rise to a 
richer set of FQH phases that extend beyond those realized in 
single-layer systems. Examples of such phases occur ?A.v = 1 
and V = 1/2. They have a rich experimental history that we 
briefly review in Sec. II. 

In this work we focus on the QHB at total Ailing factor 
V = 1/3 -f 1/3. The early experiment by Suen et al. [31] 
measured the quasiparticle excitation gap in a wide QW as 
a function of Asas- The gap was found to close around 
Asas/ ^ 0 . 1 , with an incompressible phase on either side 
of the transition. A realistic model of this system [32], that 
included EDA calculation of the band structure, reproduced 
the observed behavior of the gap. A more complete phase 
diagram as a function of both d/is and Asas/^^ was ob¬ 
tained in Ref. 33. This study, however, assumed zero width 
for each layer and was restricted to small systems. The phase 
diagram was argued to consist of three phases. Eor small djls 
and small Asas/^—> the system maintains SU(2) symmetry 
and resembles the usual ^ = 2/3 state with spin. It has been 
known that the ground state in this case is a spin-singlet ( 112 ) 
state [17, 34-36] (for an explicit wavefunction see Refs. 33 
and 37). If d/(.B is large, the layers are decoupled and the 
system is described by the Halperin (330) state, which is the 
simple bilayer Laughlin state. On the other hand, large Asas 
effectively wipes out the layer degree of freedom, and the sys¬ 
tem becomes single component. This bilayer symmetric state 
is described by the particle-hole conjugate of Laughlin’s 1/3 
wavefunction (hereafter called the 1/3 state). 

Our motivation for revisiting the problem ofi^ = 1/3-I-1/3 
QHB is twofold. Eirst, previous theoretical studies of this sys¬ 
tem have been limited to very small systems due to the expo¬ 
nential cost of exact diagonalization (ED). This limitation is 
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FIG. 1. (Color online) Phase diagram of 1/3 -F 1/3 QHB in terms of 
dimensionless layer separation d and tunnelling energy Asas- Data 
was taken with cylinder circumference L = 1Mb and layer width 
ui = 0. The dashed lines indicate sweeps performed to determine 
the nature of the phase transitions (see Sec. IV for details). Later 
in this work, additional axes will he added to this plot, driving the 
system into a non-Ahelian phase (see Fig. 7). The hlack dashed line 
and square mark the region studied experimentally in Ref. 38, and 
their observed phase transition. 


particularly severe in the present case because of the pseudo¬ 
spin degree of freedom. Recent work has demonstrated that 
to some degree this cost can be overcome by using variational 
methods such as the “infinite density-matrix renormalization 
group” (iDMRG) [39, 40]. By combining insights from ED 
and iDMRG, we are able to obtain a more accurate phase di¬ 
agram of the j/= 1/34-1/3 QHB system as a function of d 
and Asas. as shown in Eig. 1. Although our results are qual¬ 
itatively consistent with Ref. 33, the access to significantly 
larger system sizes enables us to study the order of the associ¬ 
ated phase transitions, which we And to be Arst order. 

Given that 1/3 4- 1/3 bilayer systems are experimentally 
available and allow a great deal of tunability (changing the 
layer width w, d or Asas), our second goal is to explore the 
possibility of realizing more exotic (non-Abelian) phases in 
these systems by tweaking the interaction parameters. Indeed, 
recently a number of trial non-Abelian states have been pro¬ 
posed for these systems [41-47]. At Ailing v = 1/3 4-1/3, 
the relevant candidates are the Z 4 Read-Rezayi state [48], 
the bilayer Eibonacci state [47], the “intralayer-Pfafflan” and 
“interlayer-Pfafflan” states [42]. The latter is an example of a 
spin-charge separated state and was Arst introduced in Ref. 41. 
We develop a diagnostic that detects spin-charge separation in 
the ground-state wavefunction using the entanglement spec¬ 
trum. By varying the short-range Haldane pseudopotentials in 
the bilayer system at flnite interlayer distance and tunneling, 
we And evidence for a non-Abelian phase that exhibits spin- 
charge separation and has non-trivial ground-degeneracy, con¬ 
sistent with the interlayer Pfafflan state. The phase is realized 
by either reducing the Vq or increasing the Vi pseudopotential 
component of the interaction, which may naturally occur as a 
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consequence of strong Landau level mixing. 

The remainder of this paper is organized as follows. In 
Sec. II we review some of the previous experimental work 
in QHB and QW systems. In Sec. Ill we introduce the model 
of the QHB and discuss the numerical methods and diagnos¬ 
tics for identifying the FQH phases and transitions between 
them. Sec. IV contains our main results for the phase dia¬ 
gram of 1/3 -f 1/3 QHB as a function of parameters w, d and 
Asas- We discuss in detail the three Abelian phases that oc¬ 
cur in this system, and identify the nature of the transitions 
between them. In Sec. V we explore the possible new phases 
when the interaction is varied away from the bare Coulomb 
point. We establish that the modification of short-range (Vb 
or Vi) pseudopotentials leads to a robust non-Abelian phase 
that exhibits spin-charge separation and can be identified with 
the interlayer Pfaffian state. Our conclusions are presented in 
Sec. VI. 


n. EXPERIMENTAL BACKGROUND 

In this Section we briefly review some of the important ex¬ 
periments on quantum Hall bilayers and wide quantum wells. 
As mentioned in the Introduction, one of the great advantages 
of studying these systems is the ability to experimentally tune 
parameters in the Hamiltonian, e.g., the interlayer separation 
and interlayer tunneling in a QHB. Different samples can be 
constructed with different values for these quantities. Tun¬ 
neling energy is independent of layer separation since it can 
be varied by changing the height of the potential barrier be¬ 
tween the layers without changing its width. Another conve¬ 
nient way to tune these parameters is by applying voltage bias 
to separate contacts made to each layer [49]; the variation of 
electron density p thus changes the effective fs at fixing fill¬ 
ing V via the relation p = i//27r£g. This allows d/is and 
Asas/^j^ to be tuned continuously in a single sample. 

To illustrate the typical parameter range that can be ac¬ 
cessed, we note that at = 1/2 -f 1/2 it has been possi¬ 
ble to vary djis in range 1.2^, while the interlayer tunnel¬ 
ing Asas can be either completely suppressed or as large as 
O.le^/e^B [50]. The width of individual layers in this case 
is less than d. On the other hand, in wide QWs one controls 
independently the width of the entire well and the tunneling 
amplitude Asas- The latter is defined as the energy split¬ 
ting between the lowest symmetric and antisymmetric sub¬ 
bands, and typically varies between zero and Q.2e^ jef-B- For 
systems where FQH can be observed, the physical width of 
the well is typically 30-65 nm [51]. Self-consistent numeri¬ 
cal calculations estimate that this corresponds to an effective 
bilayer distance djis = 3-7, with individual layer widths 
1.5-3 ^_b [51]. The tunability via d/ls or Asas/^j^ can en¬ 
gender new physics that does not arise in a single layer quan¬ 
tum Hall system. Two important examples of such phenomena 
have been observed to occur at total filling factors :/ = 1/2 
and V = 1. 

At total filling v = 1/2, the QHB ground state is compress¬ 
ible in the limit of both very large and very small d/ls- At 
large d/ls, it is described by two decoupled 1/4 -|-1/4 “com¬ 


posite Fermi liquids” [52] (CFL), while around d/^B = 0 it 
is the spin unpolarized 1/2 CFL. At intermediate d/lB, an 
incompressible state forms when d/iB ^ 3 [53, 54]. Numer¬ 
ical calculations performed over the years, primarily utilizing 
exact diagonalization [55-59], have confirmed that the incom¬ 
pressible state at vanishing interlayer tunneling is the Halperin 
331 state [16]. More recently, there has been some renewed 
interest in the v = 1/2 two component systems [51, 60] due 
to the possible transition into the Moore-Read Pfaffian state as 
tunneling is increased [61-63]. Analogous scenario may hold 
for QWs at filling :/ = 1/4, where the competing phases are 
the Halperin (553) state and the 1/4 Pfaffian state [64]. Very 
recently, GaAs hole systems have been shown to realize an 
incompressible state v = 1/2 near the vicinity of Landau 
level crossing [65]. 

As a second example of novel phases in QHB systems, we 
briefly mention the celebrated v = \ state (for recent reviews, 
see Refs. 50 and 66). At large d/lB the system is compress¬ 
ible (two decoupled CFLs), but undergoes a transition to an 
incompressible state for d/^B < 2, even at negligible in¬ 
terlayer tunneling. The incompressible state is represented 
by the Halperin (111) state, which can also be viewed as a 
pseudo-spin ferromagnet [67]. This wavefunction encodes the 
physics of exciton superfluidity, with an associated Goldstone 
mode [68] and vanishing of Hall resistivity in the “counter¬ 
flow” measurement setup [69, 70]. The existence of an in¬ 
compressible state (consistent with an exciton superfluid) has 
been established in numerics [71-75], though the questions 
about the details and nature of the transition, as well as the 
possibility of intermediate phases, remain open. 

The case of total Ailing v = 2/3, which is the subject of 
this paper, has been less studied compared to previous exam¬ 
ples. In the mentioned Ref. [31] the transition between a one- 
component and two-component phase was detected as a func¬ 
tion of Asas. while in Ref. [32] similar data was obtained as 
a function of the tilt angle of the magnetic held. These ex¬ 
periments have been performed on a single wide QW. More 
recently. Refs. [38] and [76] have studied v = 1/3 -f 1/3 in 
a QHB sample which directly corresponds to the model we 
study, (see Sec. Ill) By applying a voltage bias as described 
above, they perform four sweeps in the d, Asas plane. In 
one sweep [38] they And a seemingly first-order transition at 
d/(.B ~ 2, Asas/^^ ~ 0.1. This sweep, and the location of 
the observed transition, are shown in Fig. 1 . Another sweep 
entirely in the large Asas regime sees no phase transition, 
while two other sweeps are performed at small Asas- These 
sweeps see a ^ = 2/3 state at large d/fs which vanishes as 
the interlayer separation is decreased. The rest of the phase 
diagram remains to be fully mapped out. In our work we de¬ 
termine this phase diagram numerically, which can guide ex¬ 
periments towards realizing all the possible phases in this bi¬ 
layer system. Finally, we mention that very recently [77] the 
stability of fractional quantum Hall states was investigated in 
a wide quantum well system with competing Zeeman and tun¬ 
neling terms. The Zeeman splitting was controlled by an in¬ 
plane magnetic held. This system may not be fully captured 
by our model in Sec. Ill because of the potentially strong or¬ 
bital effect of an in-plane field in a wide QW. It is possible. 
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however, that the transition observed at i/ = 5/3 in Ref. 77 is 
indeed in the universality class of 1/3 (112) transition that 

we identify in Sec. IV below. 


III. MODEL AND METHOD 


A. The bilayer model 


We label the two layers of the bilayer with the index ^ £ 
{t,),}, and consider Hamiltonians of the general form 



where c''l(r) creates an electron in layer /r at the position r = 
{x, y). The first term is the Coulomb interaction, expressed in 
terms of the density operator 

n^(r) = c^l(r)c''(r). (2) 

for an electron in layer /i. The precise form of the interac¬ 
tion term depends on the details of the bilayer. The second 
term encodes tunneling between the two layers. When Vq'^ is 
SU(2) symmetric this Hamiltonian is equivalent to a z/ = 2/3 
system with spin, and in this case Asas can be thought of as 
the Zeeman splitting. 

In Eq. (1) we assumed that the perpendicular z coordi¬ 
nate has been integrated out, leading to an effective two- 
dimensional Hamiltonian. This is possible because the mag¬ 
netic field is perpendicular to the 2DEG plane, and the trans¬ 
verse component of the single body wavefunctions ip factor¬ 
izes. 


= (p:,{z±d/2)(p{Y). (3) 

The single-body wavefunctions depend on two length scales; 
the spatial separation d between the two layers in the direction 
z, and the finite layer width w of each layer. In this work we 
assume (p^ {z) is set by an infinite square well of width w. 



The Coulomb interaction in three dimensions is given by: 


V3D{x,y,z) 


is 

^yx^ + y'^ + 


(5) 


We can then recover the Coulomb interaction part of Eq. (1) 
by integrating out the perpendicular coordinate 


= j dzAz'\(p^{z)\'^\(p^{z')\'^ 

V3d[x,z- z'+ {l-5'p;)d). (6) 

Throughout this work we project the Hamiltonian (1) into the 
lowest Landau level, ignoring the effects of “Landau level 


mixing” present at finite In this case, it is possible 

to expand Vc in terms of the Haldane pseudopotentials 

which are the potentials felt by particles orbiting around one 

another in a state with relative angular momentum a. Later in 

this work we add additional Va terms to Vc in order to explore 

the neighboring phases. In experiment, such variations of the 

interaction may arise due to Landau level mixing [40, 78-85]. 

2 

Henceforth, we set the energy and length scales = 
= 1 whenever units are omitted. 


B. Numerical methods 

We work in the Landau gauge, {A^^Ay) = (y,o), 

where the single-particle orbitals with momentum = 2 ^ 
(to £ Z) are spatially localized near y = k^i^- The sys¬ 
tem is fully periodic along the x-direction, but naturally maps 
to a long-range interacting ID fermion chain along y-axis. 
We study such chains using exact diagonalization as well as 
density-matrix renormalization group [39, 40]. 

Lor the purposes of exact diagonalization (ED), it is useful 
to minimize the finite-size effects by assuming the ID chain to 
be periodic (i.e., the physical system is periodic along both x 
and y directions, or equivalently it has the topology of a torus). 
Using magnetic translation symmetry reduction of the Hilbert 
space [86], it is possible to study systems of about 10 electrons 
with pseudo-spin degree of freedom at filling 1/3 -f 1/3. The 
advantages of ED method are the direct access to the entire 
low-lying excitation spectrum, resolved ground state degener¬ 
acy, the ability to simulate complicated interactions (e.g., 3- 
body) that give rise to non-Abelian states, and compute over¬ 
laps between model wavefunctions and exact states. 

Because of the exponential cost of ED that becomes pro¬ 
hibitive for systems with pseudo-spin degree of freedom, the 
bulk of our results are obtained via the recently developed in¬ 
finite DMRG method (iDMRG) [39, 40] that allows access 
to larger system sizes. iDMRG places the Hamiltonian on an 
infinitely long cylinder of circumference L, and employs a 
variational procedure to find the ground state within the vari¬ 
ational space of matrix product states (MRS) [87-89]. MPS 
can only represent systems with a finite amount of entangle¬ 
ment S, which in turn is limited by the “bond dimension” 
X via S < log(x), while the computational resources re¬ 
quired scale as 0{x?)- In this work we used a bond dimension 
X ~ 5000-8000. On a cylinder, the entanglement scales with 
the circumference L, but is independent of the length of the 
cylinder. Therefore, while the complexity remains exponen¬ 
tial in the circumference, it is constant in the length of the 
cylinder, which provides an advantage over ED. 


C. Entanglement invariants for the identification of EQH 
phases 

All of the phases we study in this work are gapped, have 
quantized Hall conductance = |(e^//i), and have no lo¬ 
cal order parameter which can be used to distinguish between 
them. However, these phases do have different topological 
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orders, and we can therefore apply a number of recent devel¬ 
opments [39, 90-92] which demonstrate how the topological 
order of a system can be extracted from its entanglement prop¬ 
erties. 

In a topological theory, the ground state degeneracy on both 
the torus and infinitely long cylinder is equal to the number 
of anyon types. There is a special basis for the ground state 
manifold, the minimally entangled basis, in which each basis 
state I a) can be identified with an anyon type a [90, 93, 94]. 

By measuring how various entanglement properties of |a) 
scale with the circumference L, we can measure; the quan¬ 
tum dimensions da [93, 95]; the internal quantum numbers 
(spin, charge, etc.) of each anyon a; the “shift” § [96], or 
equivalently the bulk Hall viscosity [39]; the topological spins 
6a = and the chiral central charge c_ of the edge the¬ 

ory [39, 90, 92]. Below we provide a brief summary of these 
measurements in the context of FQH systems, and refer to 
Refs. 40 for a detailed discussion. 

To measure entanglement properties we divide the cylin¬ 
der in orbital space into two semi-infinite halves L/R and 
Schmidt decompose the state as jT') = \d )l ® \d)R- 

The entanglement entropy is defined as S' = — log A^. 

In ground state |a), the entropy Sa scales as [93, 95] 

Sa=/3L-log^ + 0(e-^/«), (7) 

where da is the quantum dimension of anyon a, and 27 is the 
total quantum dimension of the topological phase. The cor¬ 
rections are set by a length scale ^ which need not be directly 
related to the physical correlation length. 

To measure a U(l) charge Qa for anyon a, we partition 
the total charge operator into its components to the left / right 
of an entanglement cut, Q = Ql + Qr- The left Schmidt 
states are eigenstates of Qr, QL\p-,a)^ = Qfj,-a\p-',a)^, 
where |/i; a)^ are the Schmidt states of ground state |a) and 
G Z in units where the elementary charge is 1. The 
charge Qa of anyon a is given by the charge polarization in 
the ground state, which can be expressed as an “entanglement 
average” [39] 

g27riQ„ ^ g27riX:^ 

Qa is defined modulo 1. In the bilayer systems with U(l) x 
U( 1) symmetry we can apply the measurement for both layers 
to get two charges. 

Rotating the cylinder can also be viewed as a U(l) charge, 
whose generator is the momentum K. Its eigenvalues Ka can 
be combined with certain analytically calculable properties of 
the Landau levels to recover the Berry phase for an adiabatic 
Dehn twist (modular transformation). Similar to the charge, 
the resulting phase Ta = exp(27riMa) may be computed from 
an entanglement average: 

Ma=^ XlK^-a + analytic terms. (9) 


FQH Phase 

Ground-state 
degenacy 

S 

Spin-charge 

separation 

c 

(330) 

9 

3 


2 

( 112 ) 

3 

1 


0 

1/3 

3 

0 


0 

Z 4 Read-Rezayi [48] 

15 

3 


2 

Interlayer-Pfafflan [41] 

9 

3 

/ 

5/2 

Bonderson-Slingerland [44] 

9 

4 

/ 

5/2 

Intralayer-Pfafflan [42] 

27 

3 

/ 

3 

Bilayer Fibonacci [47] 

6 

7 


14/5 


TABLE I. Possible candidate states at u = 1/3 -I- 1/3 and their 
observed properties. We call a phase “spin-charge separated” if one 
can consistently assign charge/spin to the excitations, with one such 
excitation having neutral charge and pseudo-spin ±1/2 (see Sec. V). 


Ma is the “momentum polarization”, scaling as [39, 92] 

( 10 ) 

Here § is the shift, ha is the topological spin of anyon a, and 
c_ is the chiral central charge of the edge. 

The shift § [96] is an constant mismatch between the num¬ 
ber of flux iV$ and electrons Ne required to realize the ground 
state of the phase on the sphere, 7V$ = Ne/i' — §, and plays 
a particularly important role in our analysis. For the (330), 
(112), 1/3 states and the interlayer-Pfafflan (introduced in 
Sec. V below) the shift takes values § = 3,1,0, 3 respectively 
(see Tab. I), so distinguishes most of the phases. Because § in 
these cases is an integer and the dominant contribution to Ma, 
it converges very quickly and is far easier to measure than ha, 
c_ or da ■ 


IV. ABELIAN PHASE DIAGRAM 

In this Section we study the z/=l/3±l/3 QHB system as 
a function of experimentally relevant parameters: interlayer 
separation (d), tunelling (Asas), and layer width (w). We de¬ 
termine the phase diagram using the topological characteriza¬ 
tion explained in Sec. IIIC, and And three different Abelian 
phases [33]: decoupled = 1/3 bilayers (330) or the bilayer 
Laughlin phase, a bilayer-SU(2) symmetric spin-singlet hier¬ 
archy state (112), and a transversely polarized particle-hole 
conjugate of the Laughlin state 1 /3. 

Fig. 1 shows the phase diagram at well width w = 0 and 
cylinder circumference L = 14. Phase boundaries were deter¬ 
mined at the points marked in black; these points were found 
by performing simulations in sweeps, changing either d or 
AsaSj plotting the results. We And points where the cor¬ 
relation length and entanglement entropy have either discrete 
jumps or peaks, and we claim that these points are the phase 
transitions. The upper panels of Figs. 4, 5 and 6 show exam¬ 
ples of the correlation length data used to determine the loca¬ 
tions of these transitions. The dashed lines in Fig. 1 show the 
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sweeps where these data were taken. We note that the region 
in the vicinity of the tentative triple point is somewhat diffi¬ 
cult to resolve, but we have have not found any evidence for 
additional phases. The three Abelian phases can be intuitively 
understood in the following limiting cases. 

First, when Asas is small and d is large the two layers inter¬ 
act only weakly, and we have two decoupled Laughlin states. 
Second, when Asas is extremely large the single particle or¬ 
bitals are superpositions of both layers. Both symmetric and 
antisymmetric superpositions are possible, but when Asas is 
very large the antisymmetric superpositions are energetically 
forbidden (the energy difference between the two states is 
Asas), so we can view the system as a single quantum well 
with 1 / = 2/3, whose ground state is the particle-hole conju¬ 
gate of the Laughlin 1/3 state, which we call the 1/3 state. 
This state is particularly natural at d = 0, where the system 
is equivalent to a single layer with spin: the tunnelling term 
is a Zeeman field which spin-polarizes the system along the 
transverse direction. 

Third, when d = 0 and Asas = 0 the system is equivalent 
to a single-layer system with spin that has full SU(2) symme¬ 
try. The ground state is a (112) state [17, 33, 37]. 

The attentive reader might note that, topologically, the 
(112) and 1/3 phases are actually the same phase, in the sense 
that their K matrices are related by an SL(2, Z) transforma¬ 
tion. However, in the presence of rotational symmetry these 
phases have a different shift S, and so they are not the same 
phase. One may be concerned that in an experiment disorder 
will break the rotational symmetry and allow the (112) and 
1 /3 state to be continuously connected, but this is in fact not 
the case, as this transition has been seen experimentally both 
in wide quantum wells [31], and in single-layer systems with 
spin [97]. 


A. Determination of the phases 

We have determined the phases by using the entanglement 
invariants discussed in Sec. IIIC. First, we measure the mo¬ 
mentum polarization Ma in order to compute the shift §, 
which should take the values 3, 1 and 0 in the (330), (112) 
and 1/3 state, respectively. Fig. 2 shows the momentum po¬ 
larization at three representative points in the phase diagram. 
We plot Ma as a function of so by Eq. (10) we should get 
straight lines with a slope proportional to §. The green line 
(330) was taken at d = 1.6, Asas = 0, giving § Ri 3; the red 
line (112) was taken at d = 0.2 , Asas = 0 , giving § « 1; the 
blue line 1/3 was taken at d = 2, Asas = 0.1, giving § r; 0 . 
All of these values match those predicted for the appropriate 
phase. 

Fig. 3 shows entanglement spectra for the same points as 
those shown in Fig. 2. The counting and chirality of the 
low-lying entanglement spectra are unique to each phase, and 
as elaborated in Fig. 3 we find spectra consistent with each 
phase. 

The phase diagram in Fig. 1 was taken using an infinite 
cylinder with a circumference L = 14. To assess the finite size 



(i//167r^(’^)L^ 

FIG. 2. (Color online) Momentum polarization Ma for the represen¬ 
tative points from the phases in Fig. 1, plotted against M. 

The coefficient of proportionality is the shift S, which we can read 
off to be 3, 1 and 0 for the (330), (112) and 1/3 phase respectively, 
as expected. Data was taken at d = 1.6, Asas — O', d = 0.2, 
Asas = 0; and d = 2, Asas = 0.1 for the (330), (112) and 1/3 
phases, respectively. Values for the shift obtained from fitting the 
data are shown directly on the figure. 



Momentum Momentum Momentum 


FIG. 3. (Color online) Entanglement spectra for the phases in Fig. 1 : 
the (330) state, with counting of 1, 2, 5,... dispersing to the right; 
the 1/3 state, with counting 1,1,2,... dispersing to the left; the 
(112) state, which has an non-chiral spectra (being a convolution 
of a left and right mover). These results are in agreement with the 
predicted values for these phases. 


effects, we have measured the behavior of select cuts along 
the phase boundaries for L = 12-16. We found that the lo¬ 
cation of the (112) —>■ (330) transition changes with system 
size by d < 0.02. The (330) —1/3 and (112) —>■ 1/3 tran¬ 
sitions do move to smaller Asas at larger L, with a change 
from L : 12 16 of about 0.003. While the transition may 

continue shifting to slightly smaller Asas as L is further in¬ 
creased, at large d the change is small on the scale of the full 
phase diagram. 
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At smaller d, the critical value of Asas is fairly small at 
L = 14 and so we may be concerned that in the thermody¬ 
namic limit it is actually zero. We can test this at d = 0 
by exploiting the fact that tunneling acts as a simple Zeeman 
field in the spin realization, so the energetics can be fully de¬ 
termined by the energy difference between the (112) and 1/3 
phases at d = 0, Asas = 0. Using the additional symmetries 
at this point we can perform accurate finite-size scaling to ex¬ 
tract the energy difference in the thermodynamic limit, and 
we find that the transition occurs at Asas ~ 0.018. There¬ 
fore at least at small d, it appears that we have reached large 
enough sizes so that finite size effects do not change the loca¬ 
tion of the phase transition. Note that this system is formally 
equivalent to a z/ = 2/3 system with spin, and our value for 
the energy difference matches the numerical literature for the 
spin-polarization transition in that system. [98] 

We have also assessed the sensitivity to layer width w for 
select cuts through the phase boundary. In the upper panels 
of Figs. 5 and 6, we used dashed lines to show the correlation 
lengths at finite widths. We see that a finite layer width shifts 
the location of the (112) : (330) transition to larger d, while 
the (330) : 1 /3 transition is shifted to smaller Asas- At w = 1 
the boundaries have changed by about 10% compared to w = 
0 , so we don’t expect any qualitative differences in the phase 
diagram. 

Naturally there are many differences between the system 
we are studying numerically and those which are studied in 
experiments. In addition to the finite-size effects and our sim¬ 
plified treatment of layer width, we also neglect other fac¬ 
tors including Landau level mixing and disorder. One can 
therefore ask how relevant our data is to experiments, par¬ 
ticularly as to the quantitative locations of the phase transi¬ 
tions shown in Fig. 1. One way to address this is to compare 
to the experimental data which already exists. Ref. 38 stud¬ 
ied the (330) : 1/3 transition and found it at approximately 
d = 2, Asas = 0.1. The location of their observed transition 
is shown in Fig. 1. We obtain Asas ~ 0.07, and this gives us 
reason to believe that our data can be used as a guideline for 
future experiments. 


B. Order of the transitions 

The large system sizes accessible to our DMRG simulations 
allow us to assess the nature of the various phase transitions 
in Fig. 1. We find strong evidence that the (330) : 1/3 and 
(112) : 1/3 transitions are first order. The (330) : (112) tran¬ 
sition appears to be very weakly first order, though we cannot 
definitely rule out a continuous transition. To determine the 
order of the transition we check for discontinuities in dgJtJ, 
where g = AsaS: d tunes across the transition, as well as for 
divergences in the correlation length and discontinuities in lo¬ 
cal observables. 

The upper panel of Fig. 4 shows the (112) : 1/3 transition, 
at which the correlation length jumps discontinuously while 
remaining finite, indicating a strongly first-order transition. In 
the upper panels of Figs. 5 and 6 we show correlation lengths 



FIG. 4. (Color online) Data as a function of tunneling strength, cross¬ 
ing the (112) : 1/3 transition. The correlation length is flat except 
very close to the transition, where it is discontinuous. There is also 
a kink in the energy and in g. This is all consistent with a first-order 
transition. 

for (330) : 1/3 and (330) : (112) transitions. The correla¬ 
tion length peaks as the transition is approached, suggesting 
either a continuous or weakly first order transition. A contin¬ 
uous transition would be gapless, generating a large amount 
of entanglement which cannot be efficiently represented by 
an MPS; finite x effects then cutoff the divergent Con¬ 
sequently we would expect a strong dependence of ^ on the 
MPS bond dimension x- The different colored lines in the fig¬ 
ure correspond to increasing x, and we see that ^ increases 
with X, which could be consistent with a continuous transi¬ 
tion. However, a similar effect could be seen at a weakly first- 
order transition if x is not large enough to capture the state. 
Therefore we need other ways to determine the order of these 
transitions. 

Another approach is to look at behaviour of the energy at 
the transition point. For a first-order transition, we expect a 
kink in the energy, while for a continuous transition we expect 
the energy to vary smoothly. The middle panels of Figs. 4, 
5 and 6 show the energies near these transitions. The hrst 
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FIG. 6. (Color online) Data as a function of interlayer separation, 
crossing the (330) : (112) transition. The correlation length has a 
peak, while the energy has a kink and the g(r = 0) are jumps across 
the transition. This is indicative of a first order transition, though the 
transition is weaker compared to the others in the phase diagram. 


FIG. 5. (Color online) Data as a function of tunnelling strength, 
crossing the (330) : 1/3 transition. The correlation length has a 
peak near the transition, but this is consistent with both a first and 
second order transition. The middle panel shows the energy for both 
the (330) and 1/3 phases (see text), and as these lines are not paral¬ 
lel the systems energy has a kink. There is also a jump in g(r = 0), 
consistent with a first-order transition. 

order (112) : 1/3 transition has a clear kink in the energy. 
The (330) : 1/3 transition also appears of have a kink. The 
system also exhibits hysteresis for both the (112) : 1/3 and 
(330) : 1/3 transitions; if we initialize the system in the 1/3 
phase it will stay in that phase even if Asas is below its critical 
value. This is of course expected in a first order transition, and 
in the middle plot of Fig. 5 we plot two separate lines, which 
are the energy of the (330) and 1/3 phases (the actual energy 
of the system is whichever of these energies is lower). We 


can see that these lines are not parallel, which clearly shows 
that there is a kink in the system’s energy and therefore the 
transition is first order. At the (330) : (112) transition we 
find a very weak kink, so we tentatively conclude all three 
transitions are first order. 

It is also useful to look at the behavior of local correla¬ 
tions, such as the real space density-density correlation be¬ 
tween electrons in different layers; 

p(r) = (n'''(r)n'*'(0)) - (n'^(r)) (n'‘'(0)), (11) 

where n^{r) was defined in Eq. (2). In the (330) phase, the 
layers are uncorrelated, and this quantity should be approxi¬ 
mately zero. In the other phases, at small r the electrons repel 
and so g{r) should be negative. We can also look at the same 
correlation function in orbital space instead of real space; 

= {nl,n^} - (nl) ( 4 ), 


( 12 ) 
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For m = 0, this quantity will be negative in the 1 /3 phase, but 
it will be small in the other phases. When the above quantities 
have different values on either side of a phase transition, we 
expect them to jump discontinuously for a first-order transi¬ 
tion and to vary continuously for a second-order transition. 

We plot these quantities in the bottom panels of Figs. 4, 5 
and 6 , and see discrete jumps in all cases. Based on the re¬ 
sults of this section we can claim that all the transitions in the 
diagram are first order, with the strongest first order transition 
being the (112) : 1/3 transition. The (330) : (112) transition 
has only a slight kink in the energy and the jump in g{r) is 
smaller than the other transitions, so this is the weakest first 
order transition in the diagram. 

In Ref. 76, four experimental sweeps in our phase diagram 
were performed. Two of these sweeps had small Asas, and 
had d ~ 1.4 — 2.8. These sweeps found av = 2/3 state which 
we take to be the (330) state at large d, but below d Ki 1.8 they 
find no QH state. We believe that this is because their exper¬ 
iments were taken at layer width w/Ib ~ 2 , which would 
move the (330) : (112) transition to larger d, putting it near 
where they observe the vanishing QH state. Furthermore, we 
have found that the (330) ; (112) transition is weakly first- 
order, implying that at the transition there is a small energy 
gap. We surmise that the quantum Hall state is not observed 
in experiment because the gap is very small near the transition, 
and so the transition point is being smeared by finite tempera¬ 
ture and disorder effects. 


C. Spin polarization 

In addition to the bilayer degree of freedom electrons carry 
spin, resulting in a four-component system. Thus far we have 
assumed the spin is polarized by the external magnetic field, 
an assumption we can test with our simulations. 

The spin-polarized 1/3 phase at d = 0,w = 0 and large 
Asas is essentially a one-component system with filling 2/3, 
while the competing spin-unpolarized state is a two compo¬ 
nent (spin) system with each component having filling 1/3. 
The spin-unpolarized case has a lower Coulomb energy pro¬ 
portional to oc (this is why we find ( 112 ) in the 

equivalent bilayer problem), while the spin-polarized state 
gains a Zeeman energy proportional to the applied field B. 
For systems at fixed = 2/3, for a small perpendicular mag¬ 
netic field (and proportionally small density), the system will 
be in a spin-unpolarized state, while for large magnetic field 
(and density) the system will spin polarize. The spin base 
case been studied both numerically [98] and experimentally 
[97], but the results to not agree, with the numerics predicting 
a critical magnetic field of Ri 11 T and experiments measur¬ 
ing r; 3 T. It has been proposed that the difference between 
these values is due to the finite layer width of the samples 
[97]. We are in a position to confirm this, and indeed we find 
that increasing the layer width does decrease the critical mag¬ 
netic field, with a layer width of r^ 5 magnetic lengths being 
sufficient to bring experiment and simulation into agreement. 
Thus, in context to the bilayer set up, whether the 1/3 state is 
completely spin-polarized will depend on the bilayer separa- 
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FIG. 7. Phase diagram as a function of interlayer separation d and 
the modification of the Haldane potential 5Vb- We find that as —5Vo 
is increased, a new phase appears which we believe is a bilayer-spin 
charge separated non-Abelian phase. Data is taken with zero tunnel¬ 
ing Asas = 0 and layer width ui = 0. 


tion (d) and the strength of the magnetic field. 

For the bilayer-(112) point at d = 0, lu = 0 we compute the 
energy of an SU(4) symmetric four-component system (bi¬ 
layer + spin) with each component having filling 1/6. The 
resulting state is gapless, which means that our DMRG per¬ 
forms poorly and we can only obtain a rough estimate for the 
energy. However, it appears that the magnetic field required 
to spin-polarize the system is approximately an order of mag¬ 
nitude less than that required to polarize the 1/3 phase, so this 
phase should be spin-polarized even at small magnetic fields. 

In the large-d (330) phase, the problem reduces to decou¬ 
pled layers, and it is well known that z/ = 1/3 system spin- 
spin polarizes, so we expect this will remain true for all d into 
the ( 112 ) phase. 

Also note that experimental studies [38, 76] on this system 
have observed a spin-polarized system at all the tunnelling 
strengths and interlayer separations accessed, for magnetic 
fields Hr! 4-11 T. 


V. NON-ABELIAN PHASE 

In addition to the Abelian phases shown in Fig. 1, a num¬ 
ber of non-Abelian candidates have been proposed to appear 
in the 1/3-1- 1/3 system. These include the Z 4 Read-Rezayi 
state [48], the “interlayer-Pfaffian” (iPf) [41] and “intralayer- 
Pfaffian” states [42], and the bilayer Fibonacci state [47]. 
While we find no signature of these non-Abelian phases when 
restricting to the lowest Landau level and tuning the param¬ 
eters d, w, and AsaSj experimental samples certainly con¬ 
tain further tuning parameters we have neglected. To account 
for those, we have further perturbed the model with Haldane 
pseudopotentials Vq and Vi. Remarkably, we find that a mod¬ 
ification of the interlayer interaction, either through an attrac¬ 
tive hard core —5Vq or repulsive hollow-core SVi, is sufficient 
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FIG. 8. (Color online) Correlation length and energy for the spin- 
charge separated state as a function of SVq for d = 0.5, showing 
a clear first-order transition at SVo ~ 0.16. Note that correlation 
length increases rapidly as Vo is made further negative. 


to drive the system into a non-Abelian phase over a range of 
layer separations d. In Fig. 7 we show the phase diagram at 
hxed Asas = 0, w = 0, as we scan d and the interlayer per¬ 
turbation —SVq. We hnd that for all interlayer separations d 
it is possible to reduce Vq enough to reach a new phase con¬ 
sistent with the interlayer-Pfafhan (iPf) state, the evidence for 
which we present in this section. 

Fig. 8 shows a plot of correlation length and energy as a 
function of 5Vq for d = 0.5. There is clearly a peak in the 
correlation length and a kink in the energy at 5Vo ~ 0.16, 
indicative of a first-order phase transition. The other points 
in Fig. 7 were determined from similar data. As —SVq is in¬ 
creased much further, we see that the correlation length con¬ 
tinuously increases, and eventually the iDMRG becomes un¬ 
stable [shaded area in Fig. 7]. Based on small systems studied 
by ED, in this regime we expect a strongly-paired phase where 
electrons form tightly bound pairs in real space [14, 16]. Upon 
even further increase of —SVq [not shown in Fig. 7], using ED 
we hnd symmetry-broken, CDW and clustered phases [99]. 

In the new intermediate 6Vo phase the iDMRG hnds two 
nearly-degenerate ground states which we label |U^) and 
These states in fact triple the unit cell along the cylin¬ 
der, so by translating |n^), |fl^) we know there are at least 
six ground states in total. This must be understood as a lower 
bound on the degeneracy, as there is no general way to guar¬ 
antee iDMRG hnds all possible ground states. 

Our evidence for identifying the novel phase with the iPf is 


hve-fold. 

1. The shift is § = 3, as determined by the momentum 
polarization. 

2. Erom the ground state |0^) we deduce there is a anyonic 
excitation that carries pseudo-spin ±i yet is charge 
neutral. Hence the phase is “spin-charge separated”, 
and we call this excitation the spinon. 

3. The spinon excitation is non-Abelian, with quantum di¬ 
mension dQ 2 ft! 1.4 consistent with the iPf but not the 
intralayer-Pfafhan. 

4. The momentum polarization of the two ground states 
differ by Hqz — Iiqi ft — 0 . 21 , which corresponds to 
the difference in the topological spins of the associated 
anyons. 

5. The ground states exhibit a purely chiral entanglement 
spectra with counting that varies with charge sector. 

A summary of the possible candidates is listed in Tab. I 
These observations eliminate all other known candidates for 
the 1/3 + 1/3 system. In the following Sections, we give 
a brief description of the iPf phase (Sec. V A), compute over¬ 
laps against the model wavefunction using ED (Sec. V B), and 
present evidence for spin-charge separation (Sec. V C) and 
non-Abelian statistics (Sec. V D). 

A. The interlayer-Pfafflan state 

The iPf phase was hrst introduced and extensively dis¬ 
cussed in Ref. 41, and coined the interlayer-Pfafhan in 
Ref. 42. Similar to the Moore-Read phase relevant at iz = 5/2, 
the interlayer-Pfafhan has non-Abelian Ising anyon excita¬ 
tions, which behave like unpaired Majorana zero modes. But 
the iPf phase is even more interesting than the Moore-Read 
phase as it is “spin-charge separated”. Here we treat the two 
layers as an effective spin system and label them as f and 

The total charge is the sum Q = while the 

“pseudo-spin” is the difference — Q^). The lo¬ 

cal excitations are built up from neutral excitons and elec¬ 
trons. The neutral bilayer-excitons have Q = 0 and carry 
integral = 0,±1,±2,..., while the Q = 1 electrons 
carry = ±^. Thus local excitations obey the relation 
Q = 2S^ (mod 2), “locking” spin and charge together. In 
the iPf phase the electron can fractionalize into a neutral non- 
Abelian “spinon” carrying (5 = 0,5'^ = ^ and three non- 
Abelian “chargons” carrying Q = =0. Thus when in¬ 

cluding fractional excitations there are no constraints between 
charge and spin. 

A representative (model) wavefunction for the iPf phase is 
given by [41] 

= Pf ( ---) ^' 22 i({^},{w})- (13) 

\Xi-XjJ 

Here {z} and {ru} denote complex 2D coordinates of elec¬ 
trons in two layers, while {x} = {z,w} stands for coordinates 



























of all electrons, regardless of their layer index. The (221) state 
is defined as 

^'221 = - Zb)'^ ^(■Wa “ II(^a “ ^f>) 

a<.b a<.b a,b 

(14) 

There are nine anyon types in the iPf phase, which break 
up into three sets of three. Three of these anyons are over¬ 
all charge neutral and form the Ising theory: the trivial sector 
1 , a neutral fermion ijj which carries fermion parity but no 
charge, and the non-Abelian spinon excitation (j)s, which car¬ 
ries pseudo-spin but no charge [100]. In addition, 

threading 2tt flux quanta induces a charge Q = ^ + ^ Abelian 
anyon we denote by T*. The fusion rules are 

4>sXi’ = i’, (j)s X 4>s =+ $^ = 1. (15) 

By combining fluxes with the Ising sector, we obtain the 
nine anyon types: 



charge Q 

0 2 4 

3 3 

0 

1 (j) $2 

Co 0 

tjj 

Bh - 
2 

(j)s 4>s^ 


B. Exact-diagonalization overlaps 

Corresponding to the nine anyon types we should obtain 
nine degenerate ground states on the torus or an infinite cylin¬ 
der. Using the 3-body parent Hamiltonian [101, 102] for the 
model wavefunction in Eq. (13), we have verified this is in¬ 
deed the case on the torus. By performing exact diagonal- 
ization of this Hamiltonian, we And three ground states with 
zero momentum, each being 3-fold degenerate due to center- 
of-mass translations, which yields nine ground states in total. 
In the “thin-torus” limit [103, 104], these ground states reduce 
to the correct Tao-Thouless states, expected from Eq. (13). 

In small systems accessible by ED, the overlap with iPf 
model wavefunction becomes large in the novel phase iden¬ 
tified in Pig. 7. Por small systems up to 10 particles, we can 
obtain the complete set of exact ground states on the torus 
corresponding to Eq. 13, and overlap those with the same 
number of lowest states of the Coulomb interaction (possibly 
with some short-range pseudopotenials added). This defines 
an overlap matrix. The sum of singular values of the overlap 
matrix can serve as a rough indicator if the system is in the 
iPf phase or not. Por example, singular values close to zero 
would indicate the system being far from the iPf phase. In 
a finite system, singular values that can be considered “non¬ 
zero” are those larger than 1 /v/dimTf, where dim T-L is the 
dimension of the Hilbert space. Note that because of the in¬ 
variance under the center-of-mass translation, it is sufficient to 
restrict only to the three ground states with momentum equal 
to zero, i.e., we obtain a 3 x 3 overlap matrix. 
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FIG. 9. (Color online) Overlap between the iPf state and the ground 
state of Coulomb interaction with modified short-range pseudopoten¬ 
tials from ED. The system contains 8 electrons and 12 flux quanta, 
on a torus with a hexagonal unit cell. The color scale indicates the 
sum of singular values of the 3x3 overlap matrix defined in the 
main text, (a) The interaction is varied by changing d/is and adding 
SVq pseudopotential, (b) The interaction is varied by changing d/is 
and adding SVi pseudopotential. The same amount of 5V\ is added 
to both intralayer and interlayer Coulomb, (c) The effect of varying 
both Vo and Vi at fixed bilayer distance d — 1.5£b- Note that the 
iPf phase is located in the narrow red strip, and can be stabilized by 
either the reduction in Vo (a), the increase in Vi (b), or increase of 
both SVo and 5V\ (c). 


Pig. 9 summarizes the effect of varying short-range Vb and 
Vi components of the Coulomb interactions inferred from the 
overlap of the ground state (obtained by ED) and the model 
wavefunction, Eq. (13). We plot the sum of singular values 
of the overlap matrix between the exact ground state of the 
Coulomb interaction (with modified short-range components) 
and the iPf state. In Pig. 9(a),(b) we vary the bilayer distance 
d and add Vq (a) or Vi pseudopotential (b) to the Coulomb in¬ 
teraction. The system contains 8 electrons and 12 flux quanta 
on a torus with a hexagonal unit cell. We first note that the 
largest value of the overlap occurs in the narrow red strip, cor¬ 
responding to intermediate values of d and the reduction of 
Vb or, conversely, the increase of Vi. The non-zero overlap in 
this region suggests that the system is in the iPf phase. The 
ED result in Pig. 9(a) can be directly compared with the phase 
diagram obtained by DMRG in Pig. 7. We note that the vari¬ 
ation 5Vi in Pig. 9(b) assumes adding the same amount of 
6Vi to both intralayer and interlayer Coulomb pseudopoten¬ 
tial. Another possibility is to add SVi to interlayer Coulomb 
only. This yields a qualitatively similar result to Pig. 9(b) but 
with somewhat stronger finite-size effects. 

Pinally, in Pig. 9(c) we consider a combined effect of simul¬ 
taneously varying Vb and Vi. The starting point is Coulomb 
interaction at fixed bilayer distance d = 1.5 in the (330) phase. 
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In this case we find the iPf phase to be stabilized for positive 
SVq as well as positive 6Vi. Note that the largest overlap (i.e., 
sum of singular values of the overlap matrix) is roughly the 
same in all cases shown in Fig. 9. Although the magnitude of 
the overlap with the iPf is significant, it is relatively moderate 
(at maximum 1.8 compared to the “perfect” value of 3). The 
reason for this is the difficulty in fully resolving the complete 
set of iPf ground states in small finite systems. For example, 
finding only two out of three ground states will significantly 
reduce the overlaps in Fig. 9. This is responsible for small 
overlaps in at least part of the green region in Fig. 9, and leads 
to a somewhat narrower iPf phase compared to the DMRG 
result in Fig. 7. 
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C. Spin-charge separation 

We now demonstrate how we can extract the charges 
{Q, S^) of an anyon a from entanglement spectrum of its as¬ 
sociated ground state \a). Partition the cylinder with a cut 
along the circumference into “left” and “right” semi-infinite 
halves. Each left Schmidt state 1^8; a) of the MES |a) has 
■f /i- 

quantum numbers By coarse graining the Schmidt 

spectrum Xp-a over quantum-number sectors, we can look at 
the probability distribution Pa for charge Ql or spin to 
fluctuate to the left of the cut; 


1 = E PaiQpSl), 

(17a) 

Ql,si 

{QL)a= V PaiQL.SDQL, 

(17b) 

Ql.si 

{Sl)a= Y. Pa{QL.Sl)Sl. 

(17c) 


Ql,si 


The first equation expresses normalization. The “entangle¬ 
ment averages” in the second and third equation determine 
the charge and pseudo-spin of the anyon a (modulo local ex¬ 
citations). In Eig. 10, we have plotted this probability distri¬ 
bution in the spin-charge plane for the states |n^), Intu¬ 
itively the |n^) has a probability distribution associated with 
a completely neutral object, plus some number of electrons; 
in contrast |n^) has a probability distribution associated with 
a Q = 0, object, plus some number of electrons. 

The anyon associated with the latter ground state is what we 
identify with the spinon {(ps)- 

Our interpretation can be made rigorous by viewing the 
cylinder wavefunction as a ID fermion chain and appealing 
to the theory of ID symmetry-protected topological (SPT) 
phases. The internal symmetry group of the bilayer is G = 
(U(l) X U(l)) XI Z 2 , coming from particle conservation in 
each layer and the interchange of the two layers. As discussed 
in Sec. IIIC, any global symmetry group G can be restricted 
to the left half of the system in order to determine how it acts 
on left Schmidt states. In a ID symmetry-protected topologi¬ 
cal (SPT) phases, the symmetries G may be represented pro- 
jectively on the Schmidt states. The classification of ID-SPT 
phases is given by the distinct possible projective represen- 


FIG. 10. (Color online) Entanglement of for spin-charge separation 
in the non-Abelian phase. We plot the probability Pa{QL, Sp) for 
charges Ql,SP to fluctuate to the left of the cut in ground states 
a = Di /fl 2 . The center of this distribution gives the charge and spin 
of the anyon associated with the ground state. We see that ground 
state fli corresponds to a quasiparticle with = 0 and Q = 0, 
consistent with either the 1 or ^ sector; in the other ground state 
there is a quasiparticle with = 1/2 and Q = 0, consistent with 
the (ps sector. 


tations, which are in turned classified by the second group- 
cohomology classes [105]: 

ID G-symm. phases o elements of iT^(G, U(l)). (18) 

Eor our symmetry the group cohomology is 

i72(U(l)" xZ 2 ,U(l)) =Z 2 (19) 

meaning there are two possible SPT phases. We identify the 
two phases with “integral” (trivial) and “half-integral” (non¬ 
trivial) spin. On an open chain, the non-trivial SPT phase has 
a protected two-fold degenerate edge state that transforms pro- 
jectively under G. 

The claim is that the MES \(ps} has non-trivial ID-SPT or¬ 
der under G, while |1) and Ipj) are trivial. Now suppose we 
create a domain wall between the |1L) and |(/is) topological sec¬ 
tors. The domain wall must contain an anyon excitation cps- 
However, in the ID picture we have formed a domain wall be¬ 
tween trivial and non-trivial ID-SPT phases, which must have 
an emergent edge excitation. Eor our symmetry group G, the 
edge is two-fold degenerate, corresponding to the internal spin 
index = ±i of the spinon ps trapped there. 

Referring back to Eig. 10, we see that has a 2-fold 
degenerate probability distribution, which is a tell-tale signa¬ 
ture of a ID-SPT phase. One can explicitly check that |H^) 
is a non-trivial SPT phase under G, which is why we identify 
ps 

In summary, we have shown the state has an excitation with 
quantum numbers Q = 0,S^ = ±i, which rules out the 
(330), Z 4 Read-Rezayi, and Eibonacci phases. In light of this 
data, we find that | 0 ^) is consistent with \ps), while | 0 ^) is 
consistent with either |1) or \p). The absence of either the 
1 -family or ' 0 -family from our numerics is not terribly trou- 
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bling, as 1 and ijj have no symmetry properties which distin¬ 
guish them; even a slight energetic splitting of the topological 
degeneracy may consistently bias the iDMRG towards the lat¬ 
ter. 

The intralayer-Pfaffian also has a Q = 0, spinon, 

but it can be distinguished from the spinon of the interlayer- 
Pfaffian by its quantum dimension. 


D. Non-Abelian signatures 

In order to directly confirm the non-Abelian nature of 
the novel phase, we measure the quantum dimension of the 
spinon. In the iPf phase, the spinon has quantum dimension 
dtp^ = In contrast, the intralayer-Pfaffian phase has two 
kinds of (j)s excitations, each of which lives in only the top or 
bottom layers. In this phase the observed quasiparticle with 
spin-charge separation is a product of a spinon in each layer, 
and it therefore has quantum dimension d = 2. Our measure¬ 
ments of the quantum dimension therefore allow us to rule out 
the intralayer-Pfaffian. 

To make this measurement, we compute the difference 
in the entanglement entropy between |0^) and |0^). From 
Eq. (7) [93, 95], 

^02 (T) - (L) = Iog(dn2/do0 + (20) 

from which we obtain the ratio of quantum dimensions 
dn^/dni- Assuming |0^) corresponds to a Abelian anyon 
(dni = 1 ), and provided the finite-size effects are small 
enough, we extract the quantum dimension of the spin-charge 
separated anyon. In Fig. 11 we show the results of this sub¬ 
traction, for L = 12-17, for several different combinations 
of d, SVq and SVi. Finite-size and finite-x effects introduce 
significant systematic errors into our calculation of this quan¬ 
tity, leading to results for AS" = 5^2 — Sqi which vary from 
0.1 — 0.5 for different measurements. Though this prevents us 
from determining the quantum dimension precisely, we can 
still say that our results are consistent with (as 

shown by the blue dashed line in Fig. 11), and inconsistent 
with the intralayer Pfaffian value d = 2 (as shown by the green 
dashed line), and the abelian value of d = 1 . 

Furthermore, we can use a similar subtraction scheme to 
extract the (relative) topological spin of the spinon (17^) com¬ 
pared to the netrual (17^) via the momentum polarization 
[Eq. (10)]. Taking the difference of the momentum polariza¬ 
tions of the ground states 

Mo 2 - Mqi = hQ 2 - hQi + (mod 1), (21) 

we can extract the difference in topological spin A/i = hQ 2 — 
hQi. As shown in Fig. 11, we get Ah Ri —0.21 for a number 
of points in phase space. This is consistent with the identifi¬ 
cation 17^ = 0 s, as htf,^ ~ djjj = ^ — ^ = —0.1875 

in the iPf phase. We attribute the difference between the ob¬ 
served and expected values to finite-size and finite-x system¬ 
atic errors. (Note that /ij = 0, and thus we can conclude 
f7i 7^ 1 .) 



L/Ib 


FIG. 11. (Color online) Differences in entanglement entropy and 
momentum polarization for the two degenerate states as a function 
of circumference. Data was taken at a variety of different interlayer 
separations, SVo and 5Vi. The blue dashed lines show the expected 
values for the iPf phase, at —Spj = log v/2 and —h^j = — ^. 

Further support for our identification of ground states can 
be found in the entanglement spectrum. We first give the the¬ 
oretical orbital entanglement spectra for the ground states of 
the iPf phase, which depends on both the ground state | a) and 
the charge across the entanglement cut. (Note that so far in 
this work, we have given the entanglement spectra for only 
one value of electric charge crossing the entanglement cut, we 
chose the value of charge which has the lowest lying entangle¬ 
ment states. Henceforth we will be explicit about the charges.) 
For any of the nine MES and fixed charge {Q,S^) across the 
entanglement cut, the entanglement spectra counting follows 
one of three possible sequences. 

Si : 1,2,6,13,..., 

s^ : 1,3,8,19,..., (22) 

Sx : 1,3,8,18,.... 

For state |1L), the entanglement spectrum follows the si se¬ 
quence for even Q, and sequence for odd Q. For state 
I’ll;), the spectrum follows s^ and si for even and odd Q re¬ 
spectively. For state |0s), the entanglement spectrum always 
follows Sa- 

We attempt to match up the low-lying states of the ground 
states 117^) and |f7^) to those expected for the iPf phase, shown 
in Fig. 12. Typically one defines the “low-lying” entangle¬ 
ment states as those below the “entanglement gap”, which is 
a window devoid of states as the circumference is increased. 
In practice, at finite system size we observe multiple regions 
without states which could be called the entanglement gap, 
which makes it difficult to specify which levels should be 
counted. We have highlighted in the figures to indicate the 
presumed counting of iPf, but in full honesty other assign¬ 
ments are possible. The left panels shows the entanglement 
spectra for |f7^), with {Q,S^) = (0,0) (top) and = 

(—1, (bottom). Assuming states with — log(A^) < 5 are 
the low-lying states, we observe the counting 1,3,... and 
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FIG. 12. (Color online) Entanglement spectra for the putative iPf 
state. The left two panels show the entanglement spectra for the |fl^) 
state, for the charge sectors with the lowest-lying and second lowest- 
lying entanglement states, the counting for these states is 1,3,... and 
1, 2,..., as expected if |fl^) = |t/)). The right panel shows spectra for 
the |fl^) state. There are two degenerate charge sectors with lowest 
lying states. Here we show one example from each of the two sectors 
with the lowest-lying entanglement states, and we find counting of 
1,3,... in both, as expected if |fl^) = \4>s)- 


1,2,... for the two charge sectors, respectively. This suggests 
that |n^) = !'(/;), consistent with the momentum polarization 
data above. On the right panels, we showed the entanglement 
spectra for |0^) (for the same charges), which seem to indi¬ 
cate the counting 1,3,... regardless of charge, also consistent 
with the identification |S7^) = |^s)- 

In summary, the well-established spin-charge separation 
shown in Fig. 10 rules out all currently proposed wavefunc- 
tions besides the intralayer-Pfaffian and iPf phase. The entan¬ 
glement properties and the overlaps, though not conclusive, 
are incompatible with the intralayer-Pfaffian state, but appear 
consistent with the iPf state. 


VI. CONCLUSION 

In this work, we use iDMRG and exact diagonalization 
techniques to study a bilayer quantum Hall system with fill¬ 
ing 1/3 in each layer. We find a phase diagram in terms of 
the experimentally accessible parameters: layer separation, 
interlayer tunnelling, and layer width. We find three different 
phases: a phase with decoupled layers, a bilayer-spin singlet 
phase, and a bilayer-symmetric phase. We confirm the nature 
of these phases and study the transitions between them. 

We also explore the phase diagram for Coulomb interaction 
with modified short-range components (Vb and Vi). We find 
a non-Abelian phase over a wide region of parameter space. 
This phase has anyons which carry spin 1/2 and no charge. 
This observation, coupled with a study of additional entangle¬ 
ment properties and wavefunction overlaps, leads us to con¬ 
clude that the non-Abelian phase is interlayer-Pfaffian. Our 
data for the non-Abelian phase is inconsistent with all other 
known non-Abelian candidates. However, it is possible that a 
novel non-Abelian state could be constructed that can repro¬ 
duce our results. 

Although it is experimentally not feasible to directly mod¬ 
ify a given pseudopotential, there are many realistic ways to 
change the Coulomb interaction in a quantum Hall system, 
such as varying the chemical potential to place the system in a 
higher Landau level, introducing Landau level mixing, tilting 
the magnetic field, or screening the Coulomb potential. Since 
all of these perturbations can ultimately be expanded in terms 
of 14,’s, the phase diagram we find here may be helpful in 
guiding such studies towards a realization of the non-Abelian 
phase. 


ACKNOWLEDGMENTS 


We acknowledge helpful conversations with Eddy Ar- 
donne, Maissam Barkeshli, Parsa Bonderson, and Xiaoliang 
Qi. We are particular grateful to Jim Eisenstein, without 
whom this work would not have been started. S.G. is sup¬ 
ported by National Science Engineering Research Council 
(NSERC) of Canada. M.Z. is indebted to the David & Lu- 
cile Packard Eoundation. Z.P acknowledges support by DOE 
grant DE-SC0002140. Research at Perimeter Institute is sup¬ 
ported by the Government of Canada through Industry Canada 
and by the Province of Ontario through the Ministry of Eco¬ 
nomic Development & Innovation. R.M. acknowledges the 
Sherman Eairchild Eoundation for support. 


[1] D. C. Tsui, H. L. Stormer, and A. C. Gossard, Phys. Rev. Lett. 
48 , 1559 (1982). 

[2] R. E. Prange and S. M. Girvin, The Quantum Hall Ejfect 
(Springer-Verlag, New York, 1990). 

[3] R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 

[4] F. D. M. Haldane, Phys. Rev. Lett. 51 , 605 (1983). 

[5] B. I. Halperin, Phys. Rev. Lett. 52, 1583 (1984). 

[6] J. K. Jain, Phys. Rev. Lett. 63, 199 (1989). 


[7] J. M. Leinaas and J. Myrheim, II Nuovo Cimento B 37, 1 
(1977). 

[8] D. Arovas, J. R. Schrieffer, and F. Wilczek, Phys. Rev. Lett. 
53, 722 (1984). 

[9] R. Willett, J. P. Eisenstein, H. L. Stormer, D. C. Tsui, A. C. 
Gossard, and J. H. English, Phys. Rev. Lett. 59, 1776 (1987). 

[10] G. Moore and N. Read, Nucl. Phys. B 360, 362 (1991). 






























15 


[11] M. Greiter, X.-G. Wen, and F. Wilczek, Phys. Rev. Lett. 66, 
3205 (1991). 

[12] M. Greiter, X.-G. Wen, and F. Wilczek, Nucl. Phys. B 374 , 
567 (1992). 

[13] C. Nayak and F. Wilczek, Nucl. Phys. B 479, 529 (1996). 

[14] N. Read and D. Green, Phys. Rev. B 61 , 10267 (2000). 

[15] P. Bonderson, V. Gurarie, and C. Nayak, Phys. Rev. B 83, 
075303 (2011). 

[16] B. I. Halperin, Helv. Phys. Acta. 56, 75 (1983). 

[17] X. G. Wu, G. Dev, and J. K. Jain, Phys. Rev. Lett. 71, 153 
(1993). 

[18] X. G. Wu and J. K. Jain, Phys. Rev. B 49, 7515 (1994). 

[19] S. C. Davenport and S. H. Simon, Phys. Rev. B 85, 245303 

( 2012 ). 

[20] A. C. Balram, C. Toke, A. Wojs, and J. K. Jain, Phys. Rev. B 
91,045109 (2015). 

[21] N. C. Bishop, M. Padmanabhan, K. Vakili, Y. P. Shkolnikov, 
E. P. De Poortere, and M. Shayegan, Phys. Rev. Lett. 98, 
266404 (2007). 

[22] M. Padmanabhan, T. Gokmen, and M. Shayegan, Phys. Rev. 
B 80, 035423 (2009). 

[23] T. Gokmen, M. Padmanabhan, and M. Shayegan, Nat. Phys. 
6, 621 (2010). 

[24] K. S. Novoselov, A. K. Geim, S. V. Morosov, D. Jiang, M. 1. 
Katsnelson, 1. V. Grigorieva, S. V. Dubonos, and A. A. Firsov, 
Nature 438 , 197 (2005). 

[25] Y. Zhang, Y.-W. Tan, Ft. L. Stormer, and P. Kim, Nature 438, 
201 (2005). 

[26] X. Du, 1. Skachko, F. Duerr, A. Luican, and E. Y. Andrei, Na¬ 
ture 462 , 192 (2009). 

[27] K. 1. Bolotin, E. Ghahari, M. D. Shulman, H. L. Stormer, and 
P. Kim, Nature 462, 196 (2009). 

[28] F. Ghahari, Y. Zhao, P. Cadden-Zimansky, K. Bolotin, and 
P. Kim, Phys. Rev. Lett. 106 , 046801 (2011). 

[29] C. R. Dean, A. F. Young, P. Cadden-Zimansky, L. Wang, 
Ft. Ren, K. Watanabe, T. Taniguchi, P. Kim, J. Hone, and K. L. 
Shepard, Nat. Phys. 7, 693 (2011). 

[30] B. E. Eeldman, B. Krauss, J. H. Smet, and A. Yacoby, Science 
337, 1196 (2012). 

[31] Y. W. Suen, H. C. Manoharan, X. Ying, M. B. Santos, and 
M. Shayegan, Phys. Rev. Lett. 72, 3405 (1994). 

[32] T. S. Lay, T. Jungwirth, L. Smrcka, and M. Shayegan, Phys. 
Rev. B 56, R7092 (1997). 

[33] 1. McDonald and F. Haldane, Phys. Rev. B 53, 15845 (1996). 

[34] T. Chakraborty and F. C. Zhang, Phys. Rev. B 29, 7032 (1984). 

[35] E. H. Rezayi, Phys. Rev. B 36, 5454 (1987). 

[36] P. A. Maksym, J. Phys. Condens. Matter 1 , 6299 (1989). 

[37] J. E. Moore and F. D. M. Haldane, Phys. Rev. B 55, 7818 
(1997). 

[38] N. Kumada, D. Terasawa, Y. Shimoda, H. Azuhata, 
A. Sawada, Z. F. Ezawa, K. Muraki, T. Saku, and Y. Hirayama, 
Phys. Rev. Lett. 89, 116802 (2002). 

[39] M. P. Zaletel, R. S. K. Mong, and F. Pollmann, Phys. Rev. Lett. 
110 , 236801 (2013). 

[40] M. P. Zaletel, R. S. K. Mong, F. Pollmann, and E. H. Rezayi, 
Phys. Rev. B 91, 045115 (2015). 

[41] E. Ardonne, F. J. M. v. Lankvelt, A. W. W. Ludwig, and 
K. Schoutens, Phys. Rev. B 65, 041305 (2002). 

[42] M. Barkeshli and X.-G. Wen, Phys. Rev. B 82, 233301 (2010). 

[43] E. H. Rezayi, X.-G. Wen, and N. Read, (2010), unpublished, 
arXiv: 1007.2022. 

[44] P. Bonderson and J. K. Slingerland, Phys. Rev. B 78, 125323 
(2008). 

[45] M. Barkeshli and X.-G. Wen, Phys. Rev. B 84, 115121 (2011). 


[46] M. Barkeshli and X.-L. Qi, Phys. Rev. X 4, 041035 (2014). 

[47] A. Vaezi and M. Barkeshli, Phys. Rev. Lett. 113, 236804 
(2014). 

[48] N. Read and E. Rezayi, Phys. Rev. B 59, 8084 (1999). 

[49] J. P. Eisenstein, L. N. Pfeiffer, and K. W. West, Appl. Phys. 
Lett. 57, 2324 (1990). 

[50] J. P. Eisenstein, Ann. Rev. Cond. Mat. 5, 159 (2014). 

[51] J. Shabani, Y. Liu, M. Shayegan, L. N. Pfeiffer, K. W. West, 
and K. W. Baldwin, Phys. Rev. B 88, 245413 (2013). 

[52] B. 1. Halperin, P. A. Lee, and N. Read, Phys. Rev. B 47 , 7312 
(1993). 

[53] Y. W. Suen, L. W. Engel, M. B. Santos, M. Shayegan, and 
D. C. Tsui, Phys. Rev. Lett. 68, 1379 (1992). 

[54] J. P. Eisenstein, G. S. Boebinger, L. N. Pfeiffer, K. W. West, 
and S. He, Phys. Rev. Lett. 68, 1383 (1992). 

[55] T. Chakraborty and P. Pietilainen, Phys. Rev. Lett. 59, 2784 
(1987). 

[56] D. Yoshioka, A. H. MacDonald, and S. M. Girvin, Phys. Rev. 
B39, 1932 (1989). 

[57] S. He, S. Das Sarma, and X. C. Xie, Phys. Rev. B 47, 4394 
(1993). 

[58] V. W. Scarola and J. K. Jain, Phys. Rev. B 64, 085313 (2001). 

[59] M. R. Peterson and S. Das Sarma, Phys. Rev. B 81, 165304 

( 2010 ). 

[60] J. Shabani, T. Gokmen, Y. T. Chiu, and M. Shayegan, Phys. 
Rev. Lett. 103 , 256802 (2009). 

[61] K. Nomura and D. Yoshioka, Journal of the Physical Society 
of Japan 73, 2612(2004). 

[62] Z. Papic, M. O. Goerbig, N. Regnault, and M. V. Milovanovic, 
Phys. Rev. B 82, 075302 (2010). 

[63] M. R. Peterson, Z. Papic, and S. Das Sarma, Phys. Rev. B 82, 
235312 (2010). 

[64] Z. Papic, G. Mbller, M. V. Milovanovic, N. Regnault, and 
M. O. Goerbig, Phys. Rev. B 79, 245325 (2009). 

[65] Y. Liu, A. L. Graninger, S. Hasdemir, M. Shayegan, L. N. 
Pfeiffer, K. W. West, K. W. Baldwin, and R. Winkler, Phys. 
Rev. Lett. 112 , 046804 (2014). 

[66] J. P. Eisenstein and A. H. MacDonald, Nature 432, 691 (2004). 

[67] S. M. Girvin and A. H. Macdonald, in in Novel Quantum 
Liquids in Low-Dimensional Semiconductor Structures, S.D. 
Sarma and A. Pinczuk Eds (Wiley (New York), 1995). 

[68] 1. B. Spielman, J. P. Eisenstein, L. N. Pfeiffer, and K. W. West, 
Phys. Rev. Lett. 84 , 5808 (2000). 

[69] M. Kellogg, J. P. Eisenstein, L. N. Pfeiffer, and K. W. West, 
Phys. Rev. Lett. 93, 036801 (2004). 

[70] E. Tutuc, M. Shayegan, and D. A. Huse, Phys. Rev. Lett. 93, 
036802 (2004). 

[71] J. Schliemann, S. M. Girvin, and A. H. MacDonald, Phys. Rev. 
Lett. 86, 1849 (2001). 

[72] S. H. Simon, E. H. Rezayi, and M. V. Milovanovic, Phys. Rev. 
Lett. 91, 046803 (2003). 

[73] N. Shibata and D. Yoshioka, Journal of the Physical Society 
of Japan 75 , 043712 (2006). 

[74] G. Mbller, S. H. Simon, and E. H. Rezayi, Phys. Rev. Lett. 
101 , 176803 (2008). 

[75] Z. Papic and M. Milovanovic, Advances in Condensed Matter 
Physics 2011 , 614173 (2011). 

[76] N. Kumada, D. Terasawa, M. Morino, K. Tagashira, 
A. Sawada, Z. E. Ezawa, K. Muraki, Y. Hirayama, and T. Saku, 
Phys. Rev. B 69, 155319 (2004). 

[77] Y. Liu, S. Hasdemir, J. Shabani, M. Shayegan, L. N. Pfeif¬ 
fer, K. W. West, and K. W. Baldwin, (2015), unpublished, 
arXiv:150L06958vL 

[78] W. Bishara and C. Nayak, Phys. Rev. B 80 , 121302 (2009). 


16 


[79] A. Wojs, C. Toke, and J. K. Jain, Phys. Rev. Lett. 105, 096802 

( 2010 ). 

[80] E. H. Rezayi and S. H. Simon, Phys. Rev. Lett. 106, 116801 

( 2011 ). 

[81] Z. Papic, F. D. M. Haldane, and E. H. Rezayi, Phys. Rev. Lett. 
109, 266806 (2012). 

[82] 1. Sodemann and A. MacDonald, Phys. Rev. B 87, 245425 
(2013). 

[83] M. Peterson and C. Nayak, Phys. Rev. B 87, 245129 (2013). 

[84] S. H. Simon and E. H. Rezayi, Phys. Rev. B 87, 155426 
(2013). 

[85] K. Pakrouski, M. R. Peterson, T. Jolicoeur, V. W. 
Scarola, C. Nayak, and M. Troyer, (2014), unpublished, 
arXiv: 1411.1068. 

[86] F. D. M. Haldane, Phys. Rev. Lett. 55, 2095 (1985). 

[87] M. Fannes, B. Nachtergaele, and R. Werner, Communications 
in Mathematical Physics 144, 443 (1992). 

[88] S. Ostlund and S. Rommer, Phys. Rev. Lett. 75, 3537 (1995). 

[89] S. Rommer and S. Ostlund, Phys. Rev. B 55, 2164 (1997). 

[90] Y. Zhang, T. Grover, A. Turner, M. Oshikawa, and A. Vish- 
wanath, Phys. Rev. B 85, 235151 (2012). 

[91] L. Cincio and G. Vidal, Phys. Rev. Lett. 110, 067208 (2013). 

[92] H.-H. Tu, Y. Zhang, and X.-L. Qi, Phys. Rev. B 88, 195412 
(2013). 


[93] A. Kitaev and J. Preskill, Phys. Rev. Lett. 96, 110404 (2006). 

[94] H. Li and F. Haldane, Phys. Rev. Lett. 101, 010504 (2008). 

[95] M. Levin and X.-G. Wen, Phys. Rev. Lett. 96, 110405 (2006). 

[96] X. G. Wen and A. Zee, Phys. Rev. Lett. 69, 953 (1992). 

[97] J. P. Eisenstein, H. L. Stormer, L. N. Pfeiffer, and K. W. West, 
Phys. Rev B 41, 7910 (1990). 

[98] X. C. Xie, Y. Guo, and F. C. Zhang, Phys. Rev. B 40, 3487 
(1989). 

[99] Y. Zhang, E. H. Rezayi, and K. Yang, Phys. Rev. B 90, 165102 
(2014). 

[100] Technically the list of anyons are sets of excitations, rather 
than a single excitation. What we’ve described here are the 
properties of a representative excitation within each set, with 
all other excitations related by the addition/removal of elec¬ 
trons. 

[101] F. J. M. van Lankvelt, PhD thesis (University of Amsterdam, 
2004). 

[102] C.-H. Lee, Z. Papic, and R. Thomale, (2015), (to appear). 

[103] A. Seidel, H. Fu, D.-H. Lee, J. M. Leinaas, and J. Moore, Phys. 
Rev. Lett. 95, 266405 (2005). 

[104] E. J. Bergholtz and A. Karlhede, Phys. Rev. B 77, 155308 
(2008). 

[105] L. Fidkowski and A. Kitaev, Phys. Rev. B 83, 075103 (2011). 


